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Abstract 



A semiclassical model, based on a solution of the Vlasov equation for finite 



^! systems with moving-surface, is employed to study the isoscalar dipole modes 

' in nuclei. It is shown that, by taking into account the surface degree of 

■ freedom, it is possible to obtain an exact treatment of the centre of mass 

motion. It is also shown that a method often used to subtract the spurious 
\ strength in RPA calculations does not always give the correct result. An 

alternative analytical formula for the intrinsic strength function is derived in 
a simple confined-Fermi-gas model. In this model the intrinsic isoscalar dipole 



strength displays essentially a two-resonance structure, hence there are two 



relevant modes. The interaction between nucleons couples these two modes 
and changes the compressibility of the system, that we determine by fitting 
the monopole resonance energy. The evolution of the dipole strength profile is 
studied as a function of the compressibility of the nuclear fluid. Comparison 
with recent dipole data shows some qualitative agreement for a soft equation 
of state, but our model fails to reconcile the monopole and dipole data. 
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I. INTRODUCTION 

The isoscalar dipole mode, excited by the operator 

= J2rjY 1M (v t ), (1.1) 

i 

is particularly interesting because it can give information on the compressibility of nu- 
clear matter JD-Q] (Ref. || contains a review of experimental and theoretical studies of 
the isoscalar dipole resonance). This information is complementary to that given by the 
monopole "breathing" mode @||. 

^From the theoretical point of view the study of the isoscalar dipole mode is complicated 



by the fact that the operator ( LI ) can excite also the centre of mass (cm.) motion. If the 



model employed lacks translation invariance, the cm. excitation may appear as spurious 
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strength in the calculated response. It is well known that in a perfectly self-consistent 
random-phase approximation (RPA) the problem of spurious strength does not arise since 
the cm. strength is concentrated at zero energy, however, in practice, even a small lack of 
self-consistency can give substantial spurious dissipation at positive energy (see e.g. ||). 

Recent RPA calculations by two independent groups, using different prescriptions to 
eliminate the spurious cm. strength, give somewhat contradictory results [[<],§. The intrinsic 
isoscalar dipole strength function calculated in Ref. M typically shows two main peaks 
concentrated around 10 — 15 MeV and 20 — 30 MeV. The same feature is not shared by the 
strength functions calculated in J7|, where the low-energy peak has been almost completely 
eliminated by the subtraction of the spurious strength. More recent calculations performed 
by other groups [PlJlQfl give results that are in qualitative agreement with those of Ref. |§. 



However all RPA calculations start from a zero-order (single-particle) approximation that 
does not treat the cm. motion correctly so either perfect self-consistency is required in 
order to push all the cm. strength at zero energy, or some procedure must be employed to 
subtract the spurious strength. 

A correct description of the cm. motion is vital when evaluating the response to the 
external field ( |1.1| ), since that operator excites maily the cm. motion. Here we tackle the 
problem of separating the spurious cm. excitation from the physically interesting intrinsic 



excitation by using a semiclassical method based on the Vlasov equation. In Ref. [11 



solution of the linearized Vlasov equation for finite systems with moving surface has been 



derived, which is well suited to describe cm. motion (see also Ref. |12|]). Clearly quantum 
calculations are, in principle, more rigorous than a semiclassical one, but our hope is that 
the physical insight allowed for by the semiclassical method will make our results more 
transparent. 

We calculate the response of nuclei to an external field of the kind 

V ext (r,t) = P6(t)r 3 Y 1M (T) (1.2) 



within the model of Ref. [11|. In that model the field ( |1.2|) does not give spurious dissipation 



at positive energy because, like in a perfectly self-consistent RPA calculation, the cm. 



excitation strength is concentrated at zero energy, exactly. Moreover, in the model of ||TT 
contrary to RPA, the cm. response is treated correctly already in the starting zero-order 
approximation. 



The RPA (including its relativistic generalizations [T3-|r5|) is only one of the two main 
approaches that have been used to study the dipole compression mode, the other one being 
the fluid-dynamic approach P,|T5|. Our present kinetic-theory approach can be viewed as 



intermediate between RPA and fluid dynamics. 

II. FORMALISM 

In this paper we need to evaluate the following kind of response functions: 

K jk {u) = drr j Y 1M (r)5g k (r,u), (2.1) 

more specifically, we need the functions IZn, ^13, ^33 (we put a tilde over the moving- 
surface response functions to distinguish them from the untilded fixed-surface ones). The 
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fluctuation Sgk(r,u) is the time Fourier transform of the density fluctuation Sgk(r,t) in- 
duced by an external field V ext (r, t) = (5S{t)r k Y 1M {i). This fluctuation can be obtained by 
integrating over momentum the phase-space density fluctuation Sn k (r, p, t) that is given by 
the solution of the linearized Vlasov equation, either with fixed-surface [|17] or with moving- 
surface [|11|,|13 boundary conditions: 

Sg k (r, t) = J dp Sn k (r, p, t). (2.2) 

The solution of the linearized Vlasov equation, both with fixed and moving surface, can be 
obtained at different levels of approximation. 

Here we consider first a zero-order approximation, obtained by neglecting the nucleon- 



nucleon interaction in the bulk. Note that, while the zero-order approximation of Ref. fl7 



corresponds to nucleons moving in the static equilibrium mean field and being reflected 
by the static equilibrium surface, the zero-order approximation of [JTlJ] and |12| does take 
into account partially the density fluctuations induced by the external driving force by 
considering particle reflections from a moving surface. Thus, what has been called the zero- 
order approximation in [12] includes more dynamics than the corresponding approximation 
ofQ. 

When working with a moving surface it is convenient to introduce a modified density 
fluctuation Sg that is related to Sg in the following way. The quantities we are interested in 
are of the kind (F is an arbitrary function of coordinates) 

SF(t)=f dvF(r)g(r,t) - f drF(v)g , (2.3) 

JV(t) JVo 

where g is the equilibrium density, V the equilibrium volume, and V(t) the perturbed 
time-dependent volume. Consider the integral 

/ drF(r)g(r,t) = / drF(r)g + / drF(r)g + [ drF(r)Sg(r,t) + 0(v? xt ), 

JV(t) Jv ' JAV(t) JVo V 1 

(2.4) 

where AV(t) = V(t) — Vq; we can include the contribution of J*Av(t) drF(r)g into 
f Vo drF(r)5g(r, t) if we define 



8g(r, t) = 5g(r, t) + ^ ^(^ - R)5R(t, t), (2.5) 

with 5R(r, t) giving the fluctuation of the equilibrium radius R induced by the external 
force. Thus, when evaluating the response functions ( ]2.1|) we shall use Sg instead of Sg and 
integrate over the equilibrium volume Vq only. 

In the zero-order approximation of |12j the density fluctuation Sg(r, u) is given by the 
sum of two terms, one corresponding to the fixed-surface solution of the Vlasov equation, 
and one giving the moving-surface contribution: 

Sg°(r, u) = Sg°(r, u) + Sg° surf (r, u). (2.6) 



The structure of the last term in this equation is similar to that of the last term in Eq. fl2.5|) 
since Sg° sur j is also proportional to the fluctuation SR of the nuclear surface (cf. Eq. (3.20) 
of Ref. |l2[ ) , so it is convenient to group the two terms together and write 
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n%(uj) = n Jk (u) + s%(u) (2.7) 

for the response functions in Eq. (|2.1| ). 
Here 

Tl%{u) = J dxx 2 J dyy 2 x>y k D° L=1 (x, y, u) (2.8) 

are the zero-order fixed-surface response function calculated according to Ref. [|l7j (the 
propagator D^(x,y,uj) has been defined in [[T7|]), while the quantities 

S° k (u) = g R^5R° 1M (uj) + J dxr^Y lM {v)8e surf {v,u), (2-9) 

include the contributions both of the last term in Eq.( |2.5| ) and of the surface term $Q° SUT f in 
Eq.([2.6|), the quantity 8R® m (uj) is given in Eq.(3.21) of Ref. flI3| . According to that equation 

5Rl M {u) = -(3^\. (2.10) 

The functions xli^ ) anc ^ Xi( u; ) are given in the Appendix,where the detailed expressions of 
all the response functions needed in our calculation are also reported. 



III. CM. RESPONSE 

Before evaluating the intrinsic excitation strength associated with the external field ( [L2| 



we prove our previous statement that the model of Ref. [fLlJl , unlike that of Ref. ||17|| , gives 



a correct description of cm. motion already at zero order. From Eq. ( |A17|) we have 

Air mur 

as should be expected for a free particle of mass Am (see for example |19|). Since this 
response function has no poles for u ^ 0, it does not give spurious dissipation at positive uj. 

Another interesting response function is TZ® 3 (uj), which, in the present moving-surface 
model is given by (cf. Eq. (|A 1 8|) ) 

£° 3 H = it! 2 . (3.2) 

The function TZ\ Z {uj) has a direct physical interpretation as the displacement of the 
nuclear cm. induced by the external field ( |1.2j ): 

CM = ^aM. (3-3) 



Van Giai and Sagawa [2(| have proposed a method for removing the spurious strength 
associated with the field t 3 Yim (?). They suggested to use an effective external field of the 
kind (r 3 — t]r)YiM{r) and to determine the parameter r] by the requirement that the external 
field does not induce any displacement of the nuclear cm.: 



4 



(3.4) 



By using this prescription in a fluid-dynamic model, they determined r\ = | < r 2 >, where 
< r 2 > is the nucleus mean square radius (in our model < r 2 >= f-R 2 )- Here we do not 
need to use the effective field (r 3 — rjr) because the moving- surface model employed does not 
give any spurious dissipation, however, if we do this, we find 



(3.5) 



giving 



ft? 



11 



Using Eq. (|3.2j ) we get rj = R 2 , in agreement with [20|. However we note that this simple 
result is obtained in a model that already has the correct translation-invariance property 
and that this procedure with r] = R 2 cannot be used as a general prescription for subtracting 
the spurious strength in non-translationally-invariant models. Had we used this prescription 
within the fixed-surface model of Ref |DJ, that can be taken as an example of a case in 



which the problem of spurious dissipation appears, we would have found 

n-^M (37) 

but now the two response functions in the numerator and denominator are no longer pro- 
portional (cf. Eq. ( |A2|) ), and we would need a complex and cj-dependent parameter r\ to 
subtract the spurious strength. 

In order to allow for a more quantitative appreciation of spurious effects, in Fig.l we show 
the zero-order response function calculated in the fixed-surface approximation of Eq.flA^), 
that is with no corrections for cm. effects (dotted curve), and the same response function 
calculated for the effective operator (r 3 — R 2 r) (dashed curve). This effective operator 
drastically reduces the response (by almost 90% for A = 208), but it leaves the intrinsic 
strength at the same energy of the uncorrected response. The solid curve in Fig.l instead 
shows the result of our present moving-surface model, to be discussed more extensively 
in the next section. The quantities plotted in this and in the two following figures are 
the strength functions S(hu), related to the respective response functions by S(hu) = 
— -\mlZ{fiuj) . In the translation-invariant model the strength function includes a 5-function 
at the origin, which is not shown. The excitation energy is E = fkv. 



IV. INTRINSIC DIPOLE STRENGTH 
A. Confined Fermi gas 



After having shown that the model of Ref. [O does not introduce any spurious dissipation 



at to 7^ 0, we are now in a position to evaluate the nuclear response to the external field 
( [L.2D , without having to worry about spurious effects. The response function ftg 3 is given 
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by Eq.( |A19| ). This response function is qualitatively different from those discussed in the 
previous section. Like 1Z\ X and 7£° 3 , it does contain a part that is proportional to -j (given 
by Eq. ( |A20| ) ) and represents cm. excitation. This part does not give any strength at u ^ 0. 
In addition to the cm. response, 7Z® 3 contains also an intrinsic part given by Eq. ( |A21|) . In 



terms of the variable to the intrinsic response function ( |A21| ) reads 



3A R 4 



Art m uj 



1 



(4.1) 



with the cj-moments m^ k related to the s-moments ( A22| ) by 



m 



(4.2) 



Equation (|4.1|) does express the intrinsic response function associated with the operator ( |1 . 1|) 
in terms of response functions calculated in the underlying non-translationally-invariant 



model, the analogous expression given by the widely used prescription r 3 
real rj) is 



n%{u) 



We can easily check that, if we were allowed to make the approximation 1Z\ 
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UJ 



— Tfr) (with 



(4.3) 



which is not valid in the fixed-surface model, then Eq. ( |4.1| ) would be equivalent to Eq. ([4.3D, 
but the two formulae are different in general. The intrinsic strength function associated with 
the response function (JO) is shown by the solid curves in Figs. 1 and 2, while that given 
by Eq.( |4.3|) with rj = R 2 is given by the dashed curves in the same figures. The strength 
distributions given by the two formulae ([4.1]) and ( f4.3| )are qualitatively different, in spite of 
the fact that they give the same moments m 1 and m _1 . Actually, using the high-frequency 
expansion [BT] 



Eq. flOP gives 



2mL , n/ 1 



(4.4) 



d(%uj)Tuo 



7T 



3h 2 

IAttitl 



-AR 



(4.5) 



for the intrinsic energy- weighted sum rule fh\ ntr associated with the operator (|1.1|) , in agree- 
ment with Eq. (A.4) of Ref. |Z0| (in our model < r 4 >= |i? 4 ). 

The m~ l moment is also the same for the two response functions ( |4.1|) and (4.3). Taking 
the limit of small frequency, we get 



lim ni 



intr 



.00) 



6 AR 6 

357T KpQ 



(4.6) 



where Kfg — Q^f is the Fermi-gas incompressibility parameter. This limit is related to the 
hydrodynamic sum rule 
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m intr 



du-[- -lmK° intr (u)] 



to 



7T 



(4.7) 



by m 



(4.8) 



Having established the value of the two moments m intr and m intr , we can calculate the mean 
excitation energy as 




2m<r> 



(4.9) 



in agreement with U, where a phenomenological incompressibility parameter K was used 
instead of KpQ. Since Kpc ~ 200 MeV (for R = 1.2 Astra), while the phenomenological 
value used in || was K = 220 MeV, we see that, already at the Fermi-gas level, we have a 
resonable estimate for the mean excitation energy of the isoscalar dipole strength. What is 
unrealistic in the Fermi gas model, is the pressure P exerted by the nuclear fluid, but the 
incompressibility, defined as K = 9|p, happens to be within the range of 210 ± 30 MeV 
compatible with the monopole data |4]]. In any case we can always change the compressibility 
of the system by switching on the residual interaction between nucleons. 



B. Changing the compressibility 

So far we have treated the nucleus like a gas of non-interacting fermions confined to a 
spherical cavity with perfectly reflecting walls that are allowed to move under the action of 
the gas pressure induced by the external field, taking into account the interaction between 
nucleons changes the compressibility of this nuclear fluid. 

Rather than embarking in numerical calculations, here we give only a simple estimate of 
the main effects that are expected when the nucleon-nucleon interaction is switched on and, 
in order to keep the analytic insight allowed for by our semiclassical approach, we assume a 
very schematic, separable, effective interaction of the dipole-dipole type: 

u(r,v') =a^2rr'Y 1M (r)Y 1 * M (r'). (4.10) 

M 

The parameter !F , determining the strength of the interaction, will be related to the 
compressibility, which is the relevant parameter in this context. Then the fixed-surface 
collective response functions7^ifc(c<j), given by the solution of the RPA-type integral equation 
(5.28) of Ref. jT7| are: 



K lk {oo) = (k = 1, 3) , (4.11) 

1 — a1Zii(uj) 



and 



7 



n 33 (u) = n^) + a [n °^ (4.12) 

Now, if we replace the Fermi gas fixed-surface response functions TZ3 k (u) with the above 
collective response functions lZjh(u) in Eq. (|4.1| ), we obtain a new intrinsic response function 



3AR* 1 



1, ,2^'j 3 H i 



1 - >^ 



47rma; 2 ^ i i..2 «i>) i J" 1 J 

The energy- weighted moment rh\ ntr of this response function is still given by Eq. (|4.5|) , as 
can be easily checked by studying the large-tu limit. Thus we are confident that, even 
if the dipole-dipole residual interaction ( 4.10|) violates translation invariance, we are not 
introducing any spurious strength in our intrinsic response function (|4.13| ). 

Contrary to the first moment m 1 , the inverse moment m _1 is altered by the residual 
interaction. This can be checked by evaluating the limit 

fi A /? 6 5 

lim KintJuo) = — : n ; - [1 - — aK° n (0)1 , (4.14) 

mtrK ' 357r6e F [l -ctf^O)] 1 14 llV n ' v ; 



where 



n 2 9A R 2 . 
= --— . 4.15 

5 l07T €p 



Equation (|4.14|) can be written as 



6 AR® 

lim Kintriu) = -— r— , (4.16) 

Son K 



with 



* = & 'l -£«*?,((»■ («7) 

By following exactly the same procedure for monopole vibrations, for which we assume 
a monopole residual interaction of the kind 

u(r,r') = —r 2 r' 2 , (4.18) 
An 

we find an analogous expression for the monopole incompressibility: 

K - 6e F ^M^B u 19) 

where TZ° L=0 (0) is the zero- frequency limit of the monopole response function, analogous to 
( |4.15| ). We can determine K mon by comparison with the data on giant monopole resonances 
(see Table I) and then, by assuming K = K mon , we can study the dipole response. 
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In our simplified scheme for taking into account the residual interaction between nucleons, 
we replace the Fermi-gas intrinsic response function Q4.1Q with Eq. ( 4.1 3|) . Thus Eq.( [4.13 ) 
becomes our final expression for the intrinsic isoscalar dipole response function, it is a 
generalization of the confined-Fermi-gas expression ([4. 1|) to a confined fluid of interacting 
nucleons with incompressibility K. At first we treat K as a free parameter and study its 
effect on the dipole response function, then we determine the value of K by comparison with 
the monopole-resonance data and use this value to compare with the dipole-resonance data. 

The strength function associated with the response function ( |4.13| ) is shown in Fig. 3 
for three values of the incompressibility (K = K FG , K < K FG , and K > K FG ). When the 
system becomes softer, part of the strength shifts from the high-energy peak towards the 
low-energy peak, at about 13 MeV in Lead. The position of the two centroids is also slightly 
changed, but the main effect is the change in the relative weight of the two peaks. This 



behaviuor is similar to that found in the relativistic approach of Ref. pT3|| , the position of 
the peaks is also similar for comparable values of K. However, unlike Ref. we do not 
see any qualitative difference in the sensitivity of the two peaks to the compressibility of the 
system. 

We finally note that for K > K FG the strength moves to higher energy. 



C. Comparison with data 

In Table I we report experimental centroid energies of the giant monopole resonance 



taken from Ref. | 22|j , together with our calculated values. A reasonable agreement with 
experiment is obtained for K mon = 180 MeV. 

In Fig. 4 we compare our enegy-weighted strength functions with the recent data of Ref. 



23| for Pb, Sn and Zr. The quantity plotted in the figure is 

1 . 



7T 



y(E) = 100£ - -lmJl intr {E) )/m l intr , (4.20) 



where E = hu is the excitation energy. The solid curve shows our confined-Fermi-gas 
response, that is the energy-weighted isoscalar dipole strength for an assembly of non- 
interacting nucleons with incompressibility K w 200 MeV. Our strength integrated up 
to hu = 40 MeV exhausts about 90% of the energy- weighted sum rule fh\ ntr , while integrat- 
ing up to 100 MeV exhausts about 99% of the energy-weighted sum rule. It is important to 



note that the data overshoot this sum rule [23]. The reasons for this overshooting are not 



clear to us, since, unlike the isovector dipole, no enhancement factor of the energy-weighted 
sum rule is expected in this case pl|| . 

The experimental finding of |23j that the isoscalar dipole strength consists of two com- 



ponents is qualitatively reproduced by our semiclassical model ( at the Fermi gas level the 
two resonances correspond to two complex zeroes of the function Xi( s ) °f Eq.( A13j ) 0]). In 



this respect our results are also in qualitative agreement with the quantal RPA calculations 
of PH10| and |13|-[I5| ,that predict a two-resonance structure in the isoscalar dipole response. 



For the confined Fermi gas (solid line) the position of the two peaks is too high in energy. 
Taking into account the residual interaction (dashed line) improves the agreement with data, 
but does not completely remove the discrepancy, especially for the high-energy peak. 
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The dashed curve in Fig. 4 has been calculated for K = 180 MeV, which has been 
determined by fitting the energy of the isoscalar monopole resonance within the same model. 

Our value of K lies at the lower end of the uncertainty range indicated in Ref. Q and 
is somewhat smaller than the values reported in Ref. ||. However these values refer to 
saturated and symmetric (N = Z) nuclear matter. The small discrepancy is probably 
due to the fact that we are using a sharp-surface model, while of course experiments are 
performed on real nuclei that have a diffuse surface, thus our value of K does not refer to the 
density in the interior of nuclei, but rather to some average density that takes into account 
surface diffuseness. Also, in real nuclei N is larger than Z and, as pointed out in [Q, K 
decreases with increasing asymmetry. 

By comparing the two curves shown in Fig.4, we can say that an attractive interaction 
between nucleons, that decreases the value of K with respect to the initial value of 200 MeV, 
improves the agreement with data, but the position of the high-energy peak remains too 
high in energy, even for the rather small value K = 180 MeV, suggested by our fit to 
monopole data. Further decrease of K would shift the response towards lower energy and 
would redistribute the strength by enhancing the low-energy peak and depleting the high- 
energy one, moreover this would give too small values for the monopole-resonance centroid 
energies. The coupling of the two resonances in the isoscalar dipole response due to the 
residual interaction is a novel effect displayed by our calculation and, although the position 
of the peaks is not much affected by the compressibility, we find that the stregth associated 
with them is quite sensitive to the value of K. 

Another qualitative feature of the data that is well reproduced by our calculations is the 
A~z dependence of the peak positions [P^j . Actually, an interesting property of the response 



function (f4.13|) is that, if K/ep is independent of A, then this quantity, like the functions 
T^jk, can be expressed in terms of a universal function of the dimensionless parameter s 
defined in the Appendix: 

K int r.(s)=N 33 (A)u(s). (4.21) 

The function u(s) is the same for all spherical nuclei and its explicit expression can be easily 
derived from the equations given in the Appendix. The mass number A of the particular 
nucleus being studied enters only through the normalization factor 

9A R 6 

A/- 33 (A) = f±fL (4.22) 
1d7t ep 

and does not affect the position of the resonances in the variable s. Since s oc As, the 
scaling property ( f4.21p of the response function ( |4. results in an A" 3 dependence of the 



peak positions in excitation energy. 

In conclusion, our model is able to reproduce some of the features displayed by the data 
like the two-resonance structure of the response and the A-dependence of the peak position, 
but, even after an accurate subtraction of the spurious cm. strength for which a new formula 
has been derived, it fails to reconcile the monopole and dipole data. This failure is shared 
by our semiclassical model with the quantal RPA calculations of Ref.s [fTHlOU and |13|-15 . 
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APPENDIX: RESPONSE FUNCTIONS 



In this Appendix we give explicit expressions for the functions TCq, TZ^ 3 and TZ^ 3 . Instead 
of the variable ui, it is convenient to use the dimensionless variable 

(Al) 



v F /R ' 

where v F is the Fermi velocity, we shall use also the Fermi energy e F to characterise our 
systems. Moreover we add a vanishingly small positive imaginary part to the variable s in 
order to define the behaviour of the response functions at poles: s — > s + is. 

The fixed-surface part of the response functions ( |2.7| ) can be calculated by using the 
propagator D Q L=1 {r, r', oj) given in Ref. ||17|| . The explicit expression for A non- interacting 
nucleons at zero temperature contained in a spherical cavity of radius R is: 

n] k {s) = ^~- E E fdxx* s nN {x) QlMQ Y* ] ■ (A2) 
16vr e F n t±^ N=±1 Jo s-s nN {x) 

We have defined 

nix + Na(x) , . . 

SnN (X) = ^, A3 

X 

and a(x) = arcsm(x), while the Fourier coefficients Q* N (x) are defined according to Ref. 
17]], they are the classical limit of the radial matrix elements of the operator r k . The 



coefficients needed in our calculations are: 



QUx) = {-) nR -2~T\ (A4) 



S nN{X) 



and 



x 



2 



QIn(x) = R n tf 3 ^— -(l + ^iV^— ^ (A5) 

The functions Xk( s ) an d associated with the moving-surface part of the response 

(cf. Eq. ( |2.10|) ), are given by 



^ V- ( A-TQ k nN(x) 



X° k (s) = ~ E E / dxx* s nN (* r ' , Y ' (A6) 
8tt R 3 n ^ N = ±x h s- s nN (x) 

and 

Yi(s) = — — e F s E / dxx 2 — . (A7) 

The function is the L = 1 component of defined in Eq.(3.23) of [0, we 

have used the pole expansion of the cotangent to write it in the present form, while the 
fuctions xt( s ) gi ve the numerator of Eq. (3.21) of the same reference for external fields 
t i-y M (r). We notice the similarity between the fixed-surface response function (|A2|) and 



12 



the two functions xl( s ) an d Xi( s ) appearing in the moving-surface part of the response, in 
particular the integrand in all these functions has the same poles, while the numerators are 
different. 

By using the explicit expressions of the Fourier coefficients ([A4]) and ( |A5| ) and, repeatedly, 
the identity 

1 1 1,1 1 N 

+ — ) (A8) 



s nN s ~ s nN s s nN s ~ s nN 

it is possible to express the fixed-surface response functions in terms of the functions xl( s ) 
and A useful identity, obtained in this way, is 

|f K° lk (s) = ^(xl(s) - xg(O)) (k = 1, 3) . (A9) 

We do not report here the analogous expression for 7^33 (s) since it is too involved. Moreover 

(MO) 



3 



2e, 



3A 1 r 

IT S 2 



1 



3 
4 Jo 



1 , 9 x-^ sinfsx + Na(x)) 
dxx 2 22 



N=±l 



sin sx 



Na(x)) 



(All) 



and 



xg(o) = i? 2 x?(o) = - 



3A 

47T 



(A12) 



The advantage of expressing the fixed-surface response functions and the functions xl( s ) 
in terms of the function is due to the fact that the infinite sum over n in Xi{ s ) can 

be performed exactly, consequently we can obtain exact closed expressions for all these 
functions within this model. A useful closed expression for is 



1 3 sin(2sx) 



doc tij 



sin 2 (sx) — x 2 



(A13) 



Thus, by using the relations given above, we can connect the zero-order fixed-surface and 
moving-surface solutions of the Vlasov equation obtained in Refs. [0,0. This is particularly 
simple for and «S° 3 since 



R 4 xl(s) 



2e F s 2 

while it is somewhat more complicated for 1S33 since 



(A; = 1,3) 



(A14) 



9° 



33 V 



2e F s 2 XiO) 

By inserting xt( s ) obtained from Eq. (|A9|) into ( |A14| ), we get 



(A15) 
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3AR k+1 1 
4^ 2e F ? 



that gives 



a° (s) - ^1 



and 



(A; = 1,3) 



(A16) 



(A17) 



(A18) 



for two of the full response functions in 



For 7^33, we obtain 



^33 00 



ft 



o j 



(A19) 



with 



and 



3A R 6 1 
47r 2e F s 2 



(A20) 



K> 00 



# 6 



47r 2e F s 



i i„2 gSaW 



'13 



1 nWjjW 

2 6 £jr 



(A21) 



giving the cm. and intrinsic response, respectively. The last equation expresses the intrinsic 
response function for the operator ( |1 . 1|) in terms of non-translationally-invariant response 
functions and is the main result of the present paper. It is an exact relation within the 
present model (confined Fermi gas). The moments 



M% = / dss*[- -lvaJl%{s)\ 

JO 7T 



(A22) 



are defined in terms of the fixed-surface response functions and they can be easily evaluated 
from O. Explicitly: M\ x = §^g, M\ 3 = R 2 M\ X , and M\ 3 = ^R*M\ V 

An essential property of the intrinsic response function ( |A2l ) is that its limit for s — ► 
is finite. 
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TABLES 

TABLE I. Experimental and calculated energies of giant monopole resonance. Data from Ref. 



[22] 



Nucleus 


Exp. energy 


Calc. energy 




MeV 


MeV 


208 pk 


14.17 ±0.28 


14.1 


116 Sn 


16.07 ±0.12 


17.2 


90 Zr 


17.89 ±0.20 


18.7 
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FIGURE CAPTIONS 



Fig.l Isoscalar dipole response in fixed- and moving- surface models. The dotted line 
shows the fixed-surface response for radial dependence of operator Q(r) = r 3 , while the 
dashed line is for the effective operator Q(r) = r 3 — R 2 r. The solid line shows the zero-order 
moving-surface strength given by Eq. (|4.1| ). 



Fig. 2 The dashed and solid curves are the same as in Fig.l . Note the change of vertical 
scale. 



Fig. 3 Dipole strength for different value of incompressibility parameter. The solid curve 
shows our Fermi- gas result. 



Fig. 4 Comparison of our energy-weighted strength functions with experimental data 
of Ref. |23| for 208 Pb (a), 116 Sn (b) and 90 Zr (c). The solid curve shows the Fermi-gas 
result, corresponding to incompressibility K « 200 MeV, the dashed curve is for interacting 
nucleons with K = 180 MeV. 
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